Universal phase diagrams for the quantum spin Hall systems 
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We describe how the three-dimensional quantum spin Hall phase arises from the insulator phase 
by changing an external parameter. In 3D systems without inversion symmetry, a gapless phase 
should appear between the two phases with a bulk gap. The gapless points are monopoles and 
antimonopoles (in k space), whose topological nature is the source of this gapless phase. In gen- 
eral, when the external parameter is changed from the ordinary insulator phase, two monopole- 
antimonopole pairs are created and the system becomes gapless. The gap-closing points (monopoles 
and antimonopoles) then move in the k space as the parameter is changed further. They eventually 
annihilate in pairs, with changing partners from the pair creations, and the system opens a gap 
again, entering into the quantum spin Hall phase. 

PACS numbers: 73.43.-f, 72.25.Dc, 73.43.Nq 85.75.-d 



I. INTRODUCTION 

Spin Hall effect (SHE)^ has been attracting current 
interest, because it enables us to produce spin current 
without magnetic field or magnet. The key aspect of 
this phenomenon is that the spin current is time-reversal 
invariant, unlike the spin itself. Due to this fact, the 
spin current can be induced without breaking the time- 
reversal symmetry. The physics of spin current opens up 
a new field for the spintronics. 

In addition to the SHE in metals and doped semi- 
conductors, SHE in insulators 3 -, including the quantum 
spin Hall (QSH) effect^S, has been studied intensively. 
The quantum spin Hall phase in two dimensions (2D) is 
gapped in the bulk, while it has gapless edge modes car- 
rying spin current without breaking time-reversal sym- 
metry. The interesting point is that these edge modes 
are topologically protected. They are robust against 
weak disorder or interaction^. Although edge states 
are usually sensitive to boundary conditions such as sur- 
face roughness and impurities, the present gapless edge 
states survive even if the boundary condition is changed. 
This topological protection comes from topological order 
in the bulk, which is characterized by the Z 2 topological 
number v. v takes only two values two values v = (mod 
2) (y =even) and v = 1 (mod 2) [y =odd). v = and 
v = 1 correspond to the ordinary insulator (I) phase and 
the QSH phase, respectively. The Z 2 topological number 
v represents the number of Kramers pairs of edge states. 
This phase has been proposed theoretically in bismuth 
thin film 9 . It has also been proposed theoretically in 
CdTe/HgTe/CdTe quantum welli^, and it was demon- 
strated experimentally^. Similar effect has been pro- 
posed theoretically for three dimensions (3D) 12 i 13 , and 
is demonstrated in Bio.9Sbo.1- . The following property 
of the Z 2 topological number is important. When the 
bulk states are gapped, this topological number will not 
change as far as the interaction or nonmagnetic disorder 



is not strong enough to close the bulk gap, or to break 
to time-reversal symmetry spontaneously. 

We note that the Z 2 topological number is encoded in 
the physics of gap-closing. The Z 2 topological number is 
defined as a Pfaffian of the matrix of the time-reversal 
operator, which involves the phase of the wavefunctions 
over the whole Brillouin zone. Its calculation is involved, 
and its physical meaning is hard to understood in a intu- 
itive way. On the other hand, if we focus on the change of 
the Z 2 topological number occuring at the QSH-I phase 
transition, the change involves only the local information 
in the k space, and is much simpler. Thus by studying 
how the phase transition between the QSH and the ordi- 
nary insulating phases occur, we can get deeper insight 
into the Z 2 topological number. This transition neces- 
sarily accompanies closing of the bulk gap. We note that 
the gap closing is not so trivial as it looks. Suppose we 
change one parameter in the system and check whether 
the gap closes or not. Because of the level repulsion, in 
many cases the gap does not close due to various ma- 
trix elements for interband hybridization. These matrix 
elements should vanish simultaneously, in order to close 
the gap. In some exceptional cases the gap closes; the 
conditions for the exceptional cases are related with the 
Z 2 topological number, and these are what we pursue in 
this paper. 

In the theory of gap closing by tuning an external 
parameter, momenta which satisfy k = — k (mod G) 
play an important role, where G is a reciprocal lat- 
tice vector. Such momenta are called the time-reversal 
invariant momenta (TRIM) I\, and have the values 
r„(„ in2 „3) = (nibx + n 2 h 2 + n 3 h 3 )/2 in 3D, and 
7i={ ni n 2 ) — ( n ibi + 7i2b2)/2 in 2D, where nj — 0, 1 and 
hj are primitive reciprocal lattice vectors. The TRIM 
are the momenta with Ti = G/2 where G is a recipro- 
cal lattice vector including zero vector. The TRIM are 
crucial in the sense that they are invariant under time 
reversal. It has been revealed through the research of 
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FIG. 1: (Color online) Phase diagram for the quantum spin 
Hall (QSH) and insulator (I) phases in (a) 3D and (b) 2D. m 
is a parameter driving the phase transition, and 8 represents 
breaking of inversion symmetry. 5 = corresponds to the 
inversion-symmetric system. 

Zi topological numbers^ii^ and gap-closin g 16 i 17 that the 
situation is quite different between the systems with and 
without inversion- (Z-) symmetry 

By studying the gap-closing, one can see how the phase 
transition between the QSH and insulator phases occurs. 
We studied the 2D system in Ref. [l(| and we have ob- 
tained the universal phase diagram (Fig. [1] (b)). The 3D 
system is studied in Ref. [17| and we give a handwaving ar- 
gument to deduce the universal phase diagram in 3D, as 
shown in Fig. Q] (a). However, general theory for the 3D 
systems is still lacking, particularly for the Z-asymmetric 
systems. In this paper we describe the phase transition 
between the QSH and insulator phases in 3D and char- 
acterize its topological nature in a generic context. We 
describe the gap closing in the k space for Z-asymmetric 
systems. From the topological characterization in this 
paper, we prove that the gapless phase necessarily comes 
in between the two phases. This paper is organized as fol- 
lows. In Section |TT] we describe how the phase transition 
between the quantum spin Hall and insulator phases oc- 
curs. Section HTT1 is devoted to a calculation on the three- 
dimensional Fu-Kane-Mele model to verify the results in 
the previous section. In Section IIVI we give conclusions 
and discussions. 

We henceforth consider only clean systems without any 
impurities or disorder, while the effects of impurities and 
disorders will be discussed briefly in Section IIVI The 
time-reversal symmetry is assumed throughout the pa- 
per. Our analysis here assumes that the Hamiltonian is 
generic, and we exclude the Hamiltonians which require 
fine tuning of parameters. In other words, we exclude the 
cases which are vanishingly improbable as a real material. 

II. PHASE TRANSITION BETWEEN THE 
QUANTUM SPIN HALL AND INSULATOR 
PHASES 

As we found in Ref. Ilq , in 2D Z-symmetric systems, 
the gap closing at the QSH-I transition occurs at TRIM 
k = I\. This corresponds to the expression of the Zi 



topological number as a product of the parity eigen- 
values over all the TRIMs k = I\ over the occupied 
states^; namely at the transition the conduction and 
valence bands with opposite parities exchange their par- 
ities and the Z% topological number changes. On the 
other hand, for Z-asymmetric 2D systems, the gap closes 
at ±ko + Ti (ko ^ 0) by tuning some parameter. In 
correspondence with this gap closing, the Z2 topological 
number should be expressed as an integral over the k 
space. This is indeed the Pfaffian expression of the Z-i 
topological number— 

The phase transition in 3 D 12 i 13 can be studied sim- 
ilarly to 2D 16 . The generic phase diagram in 3D is 
shown in Fig. [lja) and is different from 2D (Fig. fjjb)). 
In Z-asymmetric 3D systems, gapless phase emerges 17 , 
which is nonexistent in 2D. This gapless phase arises 
from a topological nature of the gap-closing points 
in 3D. Namely the gap-closing points in 3D k space 
are monopoles and antimonopoles, whose "monopolc 
charges" are conserved. This conservation restricts the 
form of the QSH-I phase transition, as we see in this 
paper. 



A. General description of the QSH-I phase 
transition in 3D 

Because the QSH results from the spin-orbit coupling, 
the Zi topological number v is v = mod 2, when the 
spin-orbit coupling is zero. When we think of switch- 
ing on the spin-orbit coupling, some may undergo a 
phase transition to the QSH phase. This phase tran- 
sition changes the Z2 topological number, which means 
that it is accompanied by a closing of the bulk gap. Thus 
to search for candidate materials for the QSH phase, we 
consider tuning of a single parameter to drive the phase 
transition. Let us call the parameter m. 

As we mentioned previousl y 1 ^ 17 , the phase transi- 
tion is different whether the system considered is (i) 
Z-symmetric or (ii) Z-asymmetric. The reason for the 
difference is the following. When (i) the Z-symmetry 
is present, all the states are doubly degenerate due to 
Kramers theorem. The problem is how many parameters 
should be tuned to close the gap between the conduction 
band and the valence band, which are both doubly degen- 
erate. This number is called a codimension, and in this 
case it is five, which exceeds the number of parameters 
(k, m) . Namely there are five independent parameters for 
hybridization between the valence and conduction bands, 
and unless they are finely tuned to be zero simultane- 
ously, the gap never close, and the phase transition does 
not occur. This is interpreted as level repulsion between 
the valence and conduction bands. 

Nevertheless, there is an exceptional case in (i) Z- 
symmetric systems. At TRIM k = I\, all the states 
are classified in terms of the parity eigenvalues (= ±1). 
When the valence and conduction bands have the same 
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parities, the Hamiltonian becomes^ 

5 

ff(k) = J B (k)+^a i (k)r i , (1) 
i=l 

where a^'s and £?o are real even functions of k. Tj are 
4x4 matrices given by Ti = 1 ® t x , T 2 = a z £g> t v , 
r 3 = 1 ® r z , r 4 = cr y ® r y , and r 5 = a x ® r y , where Cj 
and Tj are Pauli matrices. The gap closes when a,(k) = 
for i = 1, • • • , 5. It means that the codimension is five, 
and the gap never closes in this case. On the other hand, 
when the valence and conduction bands have the opposite 
parities, the Hamiltonian readsil, 

4 

H (k) = Oo (k) + a 5 (k)n + £ (k)r;- , (2) 

j=l 

where eto(k) and as(k) are even functions of k, and 6") (k) 
(j = 1,2,3,4) are odd functions of k. Here T- are 4x4 
matrices given by T[ — <j z ®t x , T' 2 — l<g>T y , T' 3 — <j x ®t x , 
T' 4 = (j y ®T Xl and T' 5 = 1®t z . In this case the gap closes 
only when five equations as(k) = 0, b^{\s) = are sat- 
isfied. At the TRIM, b^\Ti) (j = 1,2,3,4) identically 
vanish, and only one condition 05 (I\) = remains to be 
satisfied. This means that the codimension is reduced 
to one. Thus if the valence and conduction bands have 
opposite parities, four matrix elements (out of five) for 
hybridization between the valence and conduction bands 
vanish identically. The resulting codimension, i.e. the 
number of parameters to be tuned for gap closing, is one. 
This is equal to the number of the tunable parameter, to. 
(We note that in this case k is fixed to be IV) To sum- 
marize the gap can close only at k = T,; for Z-symmetric 
systems. This gap-closing occurs with an exchange of 
parities between the conduction and valence bands, as is 
similar to 2D. 

On the other hand, (ii) if X-symmetry is absent, the 
bands are not degenerate (except for the points with 
k = Ti). In this case the resulting codimension D c is 
thre o 17 i 18 ' 19 , which is smaller than that in the inversion- 
symmetric case. This is because the bands are nonde- 
generate and the level repulsion is less stringent. In the 
symmetry classification of Wigner and Dyso n 20 i 21 , the 
Z-symmetry breaking makes the symmetry class from 
symplectic (D c — 5) to unitary (D c = 3). The codi- 
mension (D c = 3) is less than the number of parameters 
(m,k x ,k y ,k z ). Thus the gap can close by tuning a pa- 
rameter to. 



B. Monopole-antimonopole pair creation and 
annihilation in k space 

Henceforth we focus on (ii) the X-asymmetric systems. 
Our theory is based on the physics of gauge field in k- 
spac o 22 ' 23 . When a-th band is degenerate with another 
band at an isolated point k, such point is associated with 



a monopole for the gauge field in k-space. The gauge 
field A a (k) and the corresponding field strength B a (k) 
are defined as 

A a (k) = -i(V»a(k)|V k |Va(k)>, (3) 
B Q (k) = V k xA Q (k), (4) 

The monopole density is defined as 

p a (k) = -^V k -B Q (k) (5) 

Z7T 

Though at first sight p Q (k) vanishes identically, because 
Vk • ( Vk x ) = 0, it is not true. In some cases where the 
a-th band touches with another band at some k-point, 
the wavenumber cannot be chosen as a single continu- 
ous function for the whole Brillouin zone. In such case 
the Brillouin zone should be patched with more than one 
continuous wavefunctions2i, as is similar to the vector po- 
tential of the Dirac monopole 25 . This allows a 5-function 
singularity of p(k) at the band touching. As a result the 
monopole density has the form pQs) = qiS(k — k;) 
where qi is an integer called a monopole charge. Even 
when we vary the system by changing a parameter con- 
tinuously, the monopole charge is conserved, because it is 
quantized. The only chance for the monopole charge to 
change is to create or to annihilate a pair of a monopole 
(qi = 1) and an antimonopole (q/< — — 1). More detailed 
formulation is in Appendix A. 

In the present case, we restrict ourselves to time- 
reversal symmetric cases, where we have 

B Q (k) = -B a (-k), p Q (k) = p a (-k), (6) 

where a is the label which is a time-reversed label from 
a. It means that the monopoles distribute symmetrically 
with respect to the origin. 

At the phase transition, the gap closes between a single 
valence band and a single conduction band at k — ko . In- 
stead of considering a general Hamiltonian it is sufficient 
and much simpler to consider a 2 x 2 matrix H(k,m). 
Here we introduce an external parameter m, which con- 
trols the phase transition. We note that the following 
discussion on 2 x 2 Hamiltonian is easily generalized to an 
arbitrary Hamiltonian. The 2x2 Hamiltonian i?(k, to) is 
expanded as 

3 

H(k, m) — a (k, m) +^aj(k, m)<Ji, (7) 

i=l 

where <r, (i — 1,2,3) are the Pauli matrices. The gap 
closes when the two eigenvalues are identical, i.e. when 
the three conditions aj(k,m)=0 (i = 1,2,3) are satis- 
fied. Therefore, in general, the gap-closing point in the 
(k x , k y ,k z , TO.)-hyperspace forms a curve, which we call a 
"string". Generally, this string C occupies a finite re- 
gion in TO-direction; namely, it is vanishingly improb- 
able to lie in a single value of to. When we cut the 
string C at some value of to, the intersections in the k- 
space are the points where the gap closes, namely, the 
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monopoles and antimonopoles. Thus the string C is the 
trajectory of the monopoles and antimonopoles. Because 
the monopole charge is conserved, the monopoles and 
antimonopoles are created and annihilated only in pairs, 
which means that the trajectory C of the monopoles and 
antimonopoles forms a closed loop in the (k, to) space 
(see Fig. [5] in Appendix A) . Namely, the string C has no 
end point, because an end point of C would violate the 
conservation of monopole charge. 

Henceforth we describe how the phase transition oc- 
curs, thereby opening a gap. We consider a situation 
where one side of to, e.g. to < too is gapped while the 
other side of to, e.g. to > too is gapless. This means 
that the string C exists only in the to > too region. In 
other words, we consider an extremum of the string C . 
We pick up a gap-closing point (k, to) = (ko,rao), (i.e. 
a(ko, too) = 0) and consider the vicinity of this point. We 
investigate conditions for the point (ko,TOo) to become 
an extremum of the string C. We expand the coefficients 
to the linear order 



Oi(k, to) = ^2 MijAkj + Am, 



or in a matrix form 



a(k, to) = M Ak + AtoN, 



where Akj = kj 



koj , Am = m — mo , Mij = 



dkj 



(8) 



(9) 



and 



Ni = iSjL. If the determinant of the matrix M does not 

L urn I 

vanish, the gap-closing condition, a = (01,02,03) = 
gives 



Ak 



M _1 NAto. 



(10) 



It means that a gap-closing point moves as the parameter 
to changes, and it exists on the both sides of m = m . It 
is not the case of our interest. Therefore we conclude 



detM = det 



dai 
dk.i 



= 0, 



(11) 



which is imposed in addition to Oj = 0. Thus there are 
four conditions in total, which give a set of gap-closing 
points (ko,TOo) located at an extremum of the string C. 

We now calculate behaviors of the system in the vicin- 
ity of (ko, too). If Eq. (TTTj) holds, the matrix M has a nor- 
malized eigenvector ni with null eigenvalue: Mni = 0. 
From ni we consider two additional unit vectors n a 
(a = 2,3) to form an orthonormal basis {111,112,113}. 
We adopt this basis for the k space; 



Ak = UAp = (m 



n 2 ,n 3 ) Ap 2 



(12) 



Namely, (Api, Ap 2 , Ap 3 ) is a coordinate rotated from 
(Afci, Ak 2 , Ak 3 ). From © and (TT2]) to the linear order 
in Ak and to, we have 



where = Mrij (i = 2,3). Up to this order, the gap 
closing condition, a = 0, has no nontrivial solution in 
general, because the three vectors 112, U3, N are generally 
linearly independent. It is not the case of our interest. 
Thus we have to include the next order in Ak and Ato, 
to see whether the gap closes for Ato ^ 0; 

a = AtoN + Ap 2 u 2 + Ap 3 u 3 + ^ u^A^A^ 

i,j=l,2,3,i<j 

3 

+ UiAmApi + u(Am) 2 , (14) 

where Uy, Ui and u are vectors. We first look at the gap- 
closing point a = 0. We put Am cx A where A is small, 
and investigate the order of A for each term. As we have 
seen, if Api (i — 1, 2, 3) are of the order A, the gap-closing 
condition has no nontrivial solution. The reason is that 
the right-hand side of Eq. (|14[) has no term linear in Api . 
Hence we have to consider the quadratic term in Api , for 
which we put Api oc A 1 / 2 . Then up to O(A) we have 

AtoN + Ap 2 u 2 + Ap 3 u 3 + (Api) 2 un = 0. (15) 

The solution is given by 



/ (Api) 2 
V Ap 3 



-Ato(Q- 1 N), Q = (uu,u 2) u 3 ). (16) 



Because Ap\ > 0, the solution for this exists only when 
Am has the same sign with — (<5 _1 N)i. This means that 
on one side of Am = the system is gapped, while on 
the other side the system has gapless points, 



Ap x = ± v /-(Q- 1 N)iAm, 
Ap 2 = -(Q^ 1 N) 2 Ato, 
Ap 3 = -(Q- 1 N) 3 Ato. 



(17) 
(18) 
(19) 



Thus when Ato is changed across zero, monopole- 
antimonopole pairs are created and dissociate along the 
Api direction. The trajectory of the gapless points is as 
shown in Fig. O This is exactly the case of our pursuit: 
the point of an extremum of the gap-closing. Thus we 
have shown that for given Hamiltonian with broken I- 
symmetry, such point exists in general, and the behavior 
of the gapless points in the vicinity of this pair creation 
or annihilation is described. 

Next we consider the dispersion around the gap-closing 
point. Around the monopole the dispersion is linear in k. 
It is, however, not the case around the point of monopole- 
antimonopole pair creation. From Eq. (|14p . we can derive 
k dependence of the gap at Am = (pair creation or 
annihilation). The gap is given by E g = 2|a|. Hence, for 
Am = 0. The gap behaves as 



£ g oc Ap2,A P3 ,(Api) 



(20) 



a = Ap 2 u 2 + Ap 3 u 3 + AmN, 



(13) 



Thus the dispersion along the Api direction is quadratic, 
while that along the Ap 2 and Ap 3 directions is linear. 
The Api-direction is the direction for the monopole- 
antimonopole pair to dissociate. 
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AP 2 , AP 3 




Am>0 



FIG. 2: (Color online) Trajectory for the gapless points in 
the k space. At m — mo a monopole-antimonopole pair is 
created, and they run to the opposite directions when m is 
changed. We assume (Q _1 N)i < 0, in which gapless points 
exist only for Am > 0. 
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(b) m A 



"0 II 



k=r : 



FIG. 3: (Color online) Trajectory of the gapless points for 
(a) inversion-asymmetric and (b) symmetric systems. For (b) 
inversion-symmetric systems, the gapless point is located at 
k = Ti, and is an isolated point in the m-k space. Only at 
m — mo the system is gapless. For (a) inversion-asymmetric 
systems, on the other hand, the gapless points are created 
in monopole-antimonopole pairs at m = mi, and move in k- 
space as m is varied. The system opens a gap only by pair 
annihilation of these gapless points at m — rri2- 



C. Change in the Z2 topological numbers 

According to Ref. [ID, the Z2 topological numbers Vj 
(j = 0, 1, 2, 3) in 3D are denned as 



11213 1 



{~iy° = n ^ 

n 3 -=0,l 

( _ ir=1 , 2 ,3 = yi s ni 



(21) 

n 2 n 3 , 

(22) 



where 



y/det[w(Ti)]/Pi[w(Ti)] = ±1. 



(23) 



.-k|0|u r , 



where O is the time- 



Here w nm = 

reversal operator, and u m ^ is the periodic part of 
the Bloch wavefunction. Each phase is expressed as 
i>q; [v\V2V^], which distinguishes 16 phases. Because 
among i>i, vq is the only topological number which is 
robust against disorder, the phases are mainly classified 
by v Q . When vq is odd the phase is called as the strong 
topological insulator (STI), while if it is even it is called 
the weak topological insulator (WTI). The STI and WTI 



correspond to the QSH and I phases, respectively. The 
other indices V\, V2, and v$ are used to distinguish var- 
ious phases in the STI or WTI phases, and each phase 
can be associated with a mod 2 reciprocal lattice vector 
G Vl u 2 v 3 — ^ibi + ^2b2 + i'3b3, as was proposed in Ref. H2I. 

We can relate the shape of the trajectory ("loop") 
of the gapless points with the change in the topologi- 
cal number. To see this, we note the following. From 
Eqs. (f2Tj) and ([22)) , the Z2 topological numbers in 3D are 



defined on planes >Sj : 
Brillouin zone as 

v = YlSj, 



k • = wrii with = 0, 1 in the 



n 



r,-es 



(24) 



which are gauge invariant and have the values ±1. 

Let us take the plane for example. This affects the 
numbers Uq and Based on the theory on homotopy 
characterization of the 2D QSH phase- , one can show 
the following. An intersection of the loops with the plane 
forms a set of isolated points. They are symmetric 
with respect to k = IVi 00) 1 an d the number of points is 
even. When the number of the points is 2(2N + 1) (N: 
integer), then it accompanies the change in the Zi topo- 
logical number vq and V\. Otherwise, when the number is 
AN (N: integer), then it does not accompany the change 
in the Z2 topological numbers. To show this we note 
the result in Ref. [l3|: in 2D the Z2 topological number 
v is equal (modulo 2) to an integral of the Berry cur- 
vature inside a half of the Brillouin zone plus an extra 
term coming from "contraction" of the Brillouin zonei^. 

We can regard the slice of the 3D Brillouin zone by S{ 
as a 2D Brillouin zone , which we call D [ l ] . By this 
identification we treat the 3D Z2 topological numbers in 
the same way as in 2D. When the loop intersects the half 
of the Brillouin zone d{ 1] (within S[ 1] ) once, it means 
that at some m the monopole passes through the half of 
the Brillouin zone and changes the integral of the 

Berry curvature by unity. Thus it changes v§ and v\ as 
v$ — * v' Q = + 1 (mod 2), v\ — > v' x = v\ + 1 (mod 2). 
Therefore if 2N + 1 intersections occur within the half 



of the Z2 topological numbers vq and v\ changes 

(odd <-> even), while if 2N intersections occur, v§ and v\ 
are unchanged. This completes the proof that when a 
number of intersections between the loops and the plane 
S[ is 2(2iV + l) (N: integer), it accompanies the change 
in the Z2 topological numbers, whereas AN intersections 
involve no change in the Z2 topological numbers. 

For further investigation, we consider the following ex- 
ample. Suppose we consider simultaneous pair creations 
at to = toi, and k = ±ko + Ti. As a result we have 
two monopoles and two antimonopoles. Eventually these 
will be annihilated at k = ±kg + Tj . Then there are two 
possible cases: pair annihilation occurs (A) with chang- 
ing partners (Fig. [3Ja)) and (B) without changing the 
partners (Fig. 2]) . We can show that (A) changes the Z2 
topological numbers while (B) does not. One can see the 
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FIG. 4: (Color online) Trajectory of the gapless points for 
inversion-asymmetric systems, but without the phase transi- 
tion. 



reason in two different ways. One way to see the differ- 
ence is to consider an intersection of the "loops" with the 
planes such as . In this case, on half of the Brillouin 
zone the number of intersection points is necessarily even. 

The other way to see this difference between (A) and 
(B) is to consider switching on a perturbation which re- 
stores the I-symmetry. One may wonder whether it is 
possible to restore the I-symmetry without encounter- 
ing a phase transition. When X-symmetry-breaking term 
in the Hamiltonian is sufficiently small, it is possible, 
whereas in generic systems we cannot prove that it is 
possible. Therefore, in the following we assume that it is 
possible to restore the X-symmetry without encountering 
a phase transition. In this case, the gapless loop eventu- 
ally reduces to a point in (A) (Fig. [3Ja)). Meanwhile in 
the case (B) (Fig. 2|) the loops cannot reduce to a point 
in the X-symmetric limit. In the case (B), by restoring 
the X-symmetry, each of the two loops seems to shrink to 
a point Fi); this, however, is impossible because in X- 
symmetric case the gap closing does not occur at k ^ Ti 
due to the large codimension (= 5). Thus in (B), the 
perturbation can get rid of the gapless points completely 
from the (k, m) space, i.e. the two phases of both sides 
are identical. 

From these arguments one can see that the change of 
uq is equal to the number of loops in m-k space, modulo 
2. In the case (A) (Fig. [3Ja)), there is a single loop, and 
the change of vq is one, while in the case (B) (Fig. [4]), 
the number of loop is two, and vq is unchanged. 

In the analysis in this section, we assumed that as the 
external parameters are changed, the Hamiltonian and 
the band structure change continuously. Namely, we as- 
sumed an absence of first-order transitions. Whether or 
not first-order transitions happen depends on details of 
the system and is not solely determined by system sym- 
metries and topological order. Therefore, if first-order 
transitions are taken into account, it is no longer possible 
to discuss universal properties. For example, when first- 
order transitions are allowed, the two phases, QSH and I 
phases, can transit to each other via first-order transition 
without closing a bulk gap. This situation is realized in 
the model in Ref. [26|, where the topological QSH phase is 



realized as a phase with spontaneously broken symmetry. 

III. EXAMPLE: 3D FU-KANE-MELE MODEL 

To confirm the topological discussion in the previous 
section, we take the 3D model proposed by Fu, Kane and 
Melei^ on a diamond lattice as an example. This model 
shows a transition between STI and WTI. The model is 
written as 

H^tY.clcj+i^Xso/o 2 ) ]T cJs-(4xd?.) Ci . (25) 

Here a is the size of the cubic unit cell, t represents the 
hopping, and Aso represents the spin-orbit coupling. The 
first term represents the nearest neighbor hopping, and 
the second term is a spin-dependent hopping to the next 
nearest neighbor sites. d} 4 and d?- are the vectors for the 
two nearest neighbor bonds included in the next-nearest- 
neighbor hopping. 

This four-band model is X-reversal and time-reversal 
symmetric. It means that every eigenstate is dou- 
bly degenerate by the Kramers theorem. The doubly- 
degenerate conduction and the valence bands touch at 
the three X points, X r — (2n/a)f (r — x,y,z), and 
therefore the bulk gap vanishes. To consider the phases 
with a bulk gap, suppose one changes the nearest- 
neighbor hopping to be different for the four directions 
of nearest neighbor bonds ti (i — 1,2,3,4)^2. The sys- 
tem then opens a gap between the two doubly-degenerate 
bands, while the X-reversal and time-reversal symmetries 
are preserved. In Ref. [ID, the phase boundary is studied 
when the hopping is changed slightly from the identical 
value: = t + Sti. When we set St3 = = St^, the 
phase diagram is as shown in Fig. Ela) as a function of 
Sti and 6t2 as obtained in Ref. 12. Four phases meet at 
Sti = = St 2 , where the system becomes gapless. 

The problem of our current interest is how the phase 
boundary between the WTI and the STI changes when 
the X-symmetry is broken. In the present model, the 
simplest way to break X-symmetry is to introduce an al- 
ternating on-site energy X v , like in the 2D Kane-Mele 
(KM) model on the honeycomb lattice^. In the present 
3D case, the alternating on-site energy reduces the sys- 
tem to be similar to the zincblende structure as in GaAs. 

By introducing A„ , the symmetry of the system is low- 
ered, and an analytic calculation of the phase transition 
becomes much harder. In the present case, however, with 
a procedure explained in Appendix B, we can calculate 
how the WTI-STI phase transition changes by breaking 
the X-symmetry by the A„ term. In the X-symmetric 
(A„ = 0) case, from the phase diagram (Fig. [S^a)) we 
consider St + = St± + St2 as a parameter m driving the 
phase transition, while Sti — St 2 is fixed to be a nonzero 
value, for example, Sti — 8t 2 = O.lt. This corresponds to 
the red arrow in Fig.[5]Ja). The phase diagram in the St + - 
X v plane is as given by Fig. [5Jb). When the X-symmetry 
is broken (A„ 7^ 0), the gapless region appears in the 
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(a) (b) 




Sf, 5f + =8f,+ 5f 2 



FIG. 5: (Color online) Phase diagrams for the Fu-Kane-Mele 
model with 8tz = 0, St^ — 0. t\ and t% are the bonds along 
the [111] and [ill] directions. We put A S o = O.lt. The 
axes are in the unit of t. (a) The phase diagram in Sti-St2 
plane obtained in Ref. [TJ A„ is set as zero. Each phase is 
indexed by cubic Miller indices, following Ref. Il2l . (b) The 
phase diagram in the St+-\ v plane. A„ is newly introduced 
into the Fu-Kane-Mele model. Here 5t+ — Sti +St2, while we 
fix St- — Sti — St2 = O.lt. The arrows (red) in (a) and (b) 
correspond to the identical change in parameters. 



phase diagram. This confirms our theory in the previous 
section. As we explicitly show at the end of Appendix B, 
even when the model parameters are changed perturba- 
tively, the gapless points move but never disappear. In 
this sense the gapless phase is stable. 

To confirm our theory further, we calculate the trajec- 
tory ("string") of the gapless points in k space. As the 
parameter 5t+ is changed along the arrow in Fig. [6f b) , 
the gapless points move in k space, as shown in Fig.^a). 
As a whole, the trajectory is almost circular (but not ex- 
actly) in the k space, around the X x point. Note that 
when we gradually decrease X v , the trajectory is reduced 
to the X x point. 

The change in the Z2 topological numbers can be seen 
by counting the intersection points between the trajec- 

(n) 

tory and the planes . We choose the reciprocal lat- 
tice vectors as b x = ^(-1, 1, 1), b 2 = ^(1, -1, 1), and 
b 3 = 21(1,1,-1). The X x point (= °(b 2 + b 3 )/2 = 

Tfon)) then lies on the planes s[°\ S% , and S3 , and 
the trajectory intersects these planes twice. The trajec- 
tory does not intersect the planes SP^, S^, S^. This 
means that among the Zi topological numbers, only uq, 
U2 and ^3 changes. 

In the previous section we predicted that the dispersion 
at the pair creation and annihilation is anisotropic. To 
check this, we pick up a point of pair annihilation (Fig. [7]), 
and calculate the dispersion around this point. Indeed, 
along the directions Ni and N 2 the dispersion is linear; 
meanwhile, along the direction L which is tangential to 
the trajectory of the monopole and the antimonopole, 
the dispersion is much softer and quadratic. This agrees 
with the theory in the previous section. 




0.994< fe <1 .006 

FIG. 6: (Color online) (a) Trajectory of the gapless points 
in k space, as we change the parameter St+ with X v fixed. 
The wavenumber k is shown in the unit of (2n/a). The solid 
and broken curves are the trajectories for the monopoles and 
antimonopoles, respectively. This corresponds to the arrow 
in the phase diagram in (b). 




-0.04 



FIG. 7: (Color online) (a) Dispersion around the gapless 
points when the monopole-antimonopole pairs annihilate, for 
the Fu-Kane-Mele model. The direction L is along the trajec- 
tory of the monopole and the antimonopole, while Ni and N2 
are perpendicular to L. These directions are shown in (b) in 
the trajectory for the gapless point in the k space. Energy is 
shown in the unit t, and the wavenumber is in the unit 2n/a. 



IV. CONCLUSIONS AND DISCUSSIONS 

In this paper, we described the generic phase diagrams 
involving the quantum spin Hall and insulator phases, 
by tuning an external parameter in 2D and 3D. In X- 
asymmetric 3D systems, there lies a finite region of the 
gapless phase in the phase diagram. This was checked 
for the Fu-Kane-Mele model. We described the phase 
transition in terms of the motion of the gap-closing points 
(i.e. monopoles and antimonopoles) in k space. The 
gapless phase in the X-asymmctric 3D system originates 
from the conservation of "monopole charge" 

We also studied the dispersion around the gapless 
points. In general the energy dispersion is linear around 
the gapless points. Meanwhile, around the gapless points 
at monopole-antimonopole pair creation/annihilation, 
the dispersion becomes quadratic in one direction of k, 
which is tangential to the monopole trajectory. 

It is interesting to compare the physics of gap-closing in 
the QSH phase and that in the quantum Hall (QH) phase. 
In the QH phase the bulk is gapped, and is characterized 
by the Chern number, by which this phase is distinct from 
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an ordinary insulator. For the QH phase, the behavior of 
monopoles and antimonopoles are quite similar to that 
presented so far in this paper. This point can be illus- 
trated by using a 2D model of the quantum (anomalous) 
Hall effect proposed by Haldane^l. The phase diagram is 
shown in Ref. 1271 . where the transition between the QSH 
and the insulator phases occurs by changing the model 
parameter M. This phase transition occurs when the 
gap closes only at one wavenumber in the Brillouin zone, 
which is similar to the case of the 2D QSH phased In 
3D, the gapless phase is expected to occur between the 
QH phase and the ordinary insulator phase (and also be- 
tween two QH phases with different Chern numbers). It 
is because the symmetry class is unitary, and it is simi- 
lar to the QSH phase without Z-symmetry. The 3D QH 
phase is characterized by three Chern numbers^. There- 
fore, for the transition between phases with different sets 
of Chern numbers, a gapless phase is expected to ap- 
pear in between. This gapless phase is described by a 
gapless loop C in the (to, k) space, as in the 3D QSH 
case. The topology of the loop C relative to the crystal- 
lographic directions determines the change in the three 
Chern numbers at the transition. By comparing the cases 
of the QSH phase and the QH phase, they are different 
in the following aspects. First, the QH effect is with- 
out time-reversal symmetry. This makes the physics at 
k and — k independent, in contrast with the QSH case. 
Therefore, the presence or absence of the Z-symmetry is 
inessential in the QH system. Second, because the QH 
system usually requires a strong magnetic field, and the 
motion parallel to the magnetic field is usually gapless, 
the QH phase in 3D is not easily realized in real sys- 
tems. Meanwhile there is no such constraint for the QSH 
phase, and 3D QSH phase is easily realized. Thus the 3D 
gapless phase due to the topological nature of monopoles 
is more realistic in the QSH phase than in the quantum 
Hall phase. 

When disorder is introduced in the system, the phase 
transition in three dimensions will have a gapless re- 
gion as a function of the external parameter to. In 3D, 
only for Z-symmetric systems, the disorder effect on the 
phase transition between the QSH and I phases has been 
studied^, while the Z-asymmetric systems are left to be 
analyzed. By regarding the whole system as one unit 
cell ( "supercell" ) , the similar discussion holds. We im- 
pose the periodic boundary conditions, but with allowing 
additional phase twisting for the individual directions as 
9 X , By and Z . This phase twisting plays the role of the 
wavenumbers k x , k y and k z . In this case, as the disor- 
der generally breaks the Z-symmetry, the disorder brings 
about a gapless phase. This is an interesting question, 
and is beyond the scope of the present paper. 
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APPENDIX A: GENERAL DESCRIPTION OF 
THE GAP-CLOSING POINTS IN m-k SPACE 

In this Appendix we assume that the bands are non- 
degenerate almost everywhere in k space. It allows exis- 
tence of isolated k points with band degeneracy; mean- 
while, the degeneracy in an extended region in k space, 
such as Kramers degeneracy in systems with Z- and 
time-reversal symmetry is excluded. The vectors A Q (k) , 
B Q (k) are defined in the three-dimensional k space. To 
study behaviors of the gap-closing point by changing to, 
it is convenient to consider a four-dimensional (to, k) 
space. Let us write 

fco = rn, (Al) 
and we define the following 4-vectors 

Q 

A a .i(k) ee A a .i(m,k) = -i(^ a (m, k)| — |^ a (m,k)), 

(A2) 

d d 
B a ,ij(k) ee B a ij (to, k) = —Aaj(m,k) - —A aii (m,k) 



dk; 



(A3) 



where k — (ko,ki,k2,ks) = (to, k) and i — 0,1,2,3. 
This corresponds to the vectors in Eqs. ©H]) as 
A Q (k) = (A 1 ,A 2 ,A 3 ), B Q (k) = {B 23 ,B 3U B 12 ). We 
omit the band index a henceforth, unless necessary. 
The monopole density p(k) becomes a 4-vector p — 
(po,pi,p2,Pa), where 

h'^w.ir^- (A4) 



and eujk is the totally antisymmetric tensor with 60123 
1 . We note that p(k) = p . From Eq. (|A4|) . if the band 
a is not degenerate with other bands at (m, k), Ai(k) 
is analytic, and pi(m,\<.) identically vanishes. This does 
not apply when the band is degenerate with other bands 
at some k points; at such degenerate points, p\ has a S- 
function singularity, as we see below. Such points form a 
curve in (to, k) space because the gapless condition con- 
sists of three equations, i.e. the codimension is thre o 18 i 19 
in this case. We call this gapless curve as a "string". 
We will see that the 4-vector p describes a current inside 
the string, and the total "current" inside the string is an 
integer. In general this "current" is unity. 

To see how the singularity appears, we show the fol- 
lowing; (a) a surface integral of pi over a closed three- 
dimensional hypersurface V in (to, k) space is zero (i.e. 
p is divergence-free) , and (b) a surface integral of pi over 
an open 3-dimensional hypersurface V in (to, k) space is 
quantized. To see (a) we use the Gauss theorem. 



da l3k 



eiijkpl 



1 

2^ 



i)V 



dajk w Ak = 0, (A5) 
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because dV is null. This proof relies on the fact that 
Bjk — af~^fe ~ af~^i ^ s § au S e invariant. Next we show 
(b). 



1 
3! 



da l3k eu jk pi 



1 

2^ 



da jk 



d 



A, 



(A6) 



This quantity is an integer, representing a monopole 
charge inside dV. To show this we use the Stokes theo- 
rem to this expression. In general, when the closed two- 
dimensional surface dV encloses a degeneracy point, the 
field Ak cannot be expressed as a single function on the 
surface dV . Thus the Stokes theorem can only be applied 
after dividing the surface dV into pieces, on each of which 
the wavefunction (and the field Ai) is continuous 2 ^. The 
resulting formula is expressed in terms of a phase differ- 
ence between the neighboring pieces. Because the phase 
of the wavefunctions allows a difference of a multiple of 
27r, this quantity becomes an integer. Thus from (a)(b), 
we have shown that pi describes a (divergence-free) "cur- 
rent" in the string, with its current being quantized in- 
side the string. This kind of discussion with patching of 
the wavefunction is discussed in the context of the Berry 
phas o 22 ' 23 , and also in the context of the quantum Hall 
effect, and magnetic superconductor—. 
The expression of pi is therefore given by 



Pi = 



ds dKi{s) d(4){k _ K{s)) (A7) 
as 



where K(s) = (K Q (s), K^s), K 2 (s), K 3 (s)) = 
(M(s),K(s)) describes a trajectory (i.e. string) 
of the gap-closing point in the 4-dimensional 
space, s is a parameter along the string, and 
5^{k - K(s)) = lfi=a S ( k i ~ K i( s ))- The divergence 
becomes zero because 



dpi 



ds^-5 (i) (k~ K(s)) = 0. 
as 



(A8) 



The 0-th component of (|A7p gives 
p(m, k) = J ds^^-S{m - M(s))6 (3) (k - K(s)) 
= E s S n (^) * (3) (k-K( Sl )), (A9) 



Si\ra—M{si) 

where the summation is taken over s = St which satisfies 
to = M (s^ for given to. By equating this with 

p(m,k)= 9^ (3) (k-K( Sj )), (A10) 

Si:ra—M{si) 

we get the monopole charge to be qi = sgn ^ ^ 
This means that when the "current" in the string is going 



in the increasing direction of to, it appears as a monopole 
(q = 1), whereas the decreasing direction of to corre- 
sponds to an antimonopole (q — —1), as shown in Fig. [8] 



m 




monopole 



antimonopole ', monopole 









pair 




pair 




pair 


creation 




annihilation 




creation 



FIG. 8: The string which is a set of gapless points in the 
(m, k) space. The arrow describes the direction of the flow 
vector p. 

We now return to the present case of the transition be- 
tween the QSH and the insulator phases. Recall that the 
string is confined in a restricted region toi < to < to-2 
in the to direction. Therefore the "string" becomes a 
"loop". Then it follows that at to = toi the "string" 
changes its direction in the to space, which appears as a 
pair creation of monopole and antimonopole. Similarly 
at to = to 2 a pair annihilation of monopole and anti- 
monopole results. 

In the particular cases with time-reversal invariance, 
as we are interested in, the pair creation and pair anni- 
hilation appears symmetrically with respect to k = Tj. 
Therefore, we expect that the number of loops should be 
one in the simplest case will be as shown in Fig. [3J and 
it is indeed realized in the Fu-Kane-Mele model with the 
A„ term. 



APPENDIX B: CALCULATION OF THE 
GAPLESS PHASE IN THE FU-KANE-MELE 
MODEL WITH STAGGERED ON-SITE 
POTENTIAL 

The Hamiltonian matrix for the the Fu-Kane-Mele 
model with staggered on-site potential is written as 

V F*l - Ei=i -fto 

(Bl) 

where ai are the Pauli matrices, and 1 is the 2x2 identity 
matrix. The coefficients are given by 

F = tie ia{k » +k * )/2 + t 2 e ia{k * +k * )/2 + t 3 e ia( - k * +k y )/2 + t 4 , 

(B2) 

k z a 
2 2 

and F y , F z are given similarly to F x after cyclic per- 
mutation of the subscripts x, y and z. We write F = 



,. , . kya 
b x = — 4Aso sin — — cos — cos — r— 



(B3) 
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(F x , Fy, F z ) for brevity. We study when and how the 
gap closes as the parameters change. Because the gap is 
between the second and the third bands, we need a con- 
dition when the second and third bands have identical 
eigenenergies. As the codimension is three in this case, 
the condition should be expressed as three equations, de- 
termining a string in (m, k) space. First we note that the 
spectrum of H is symmetric with respect to E = 0, as 
follows from 



(B4) 



where the superscript * denotes matrix transposition, and 
the unitary matrix U is given by 



U = a y <g> 1 = 



/ 

i 

\ i 



(B5) 



Thus the gap closes if and only if one of the eigenvalues 
of H vanishes, namely, DetH(k) = 0. This renders to 

|F | 4 + 2|F | 2 (A 2 + |F| 2 ) + (A 2 - |F| 2 ) 2 = 0, (B6) 

namely, 

ReF = 0, ImFo = 0, A 2 = |F| 2 . (B7) 



For given model parameters ti, A„, the coupled equations 
(|B7j) gives gapless points in k-space, if any. Such gapless 
points are stable against small changes of parameters, 
as we can see as follows. For brevity, let us write the 
three conditions in Eq. (|B7|) as #j(k,t) = (i = 1,2,3), 
where t denotes the set of model parameters. Let ko be 
the wavenumber of one of the gapless points. When the 
model parameters are changed, the values of gi change 
accordingly: gi — > gi + Sgi . The gapless points are then 
expected to move (k = ko — > k = ko + (5ko). Because the 
gapless conditions are expanded as 



V k g l (k ) • Sk Q + Sgi(ko) = 0, 



(B8) 



there always exists i5ko satisfying these three coupled lin- 
ear equations. Thus a perturbative change in system pa- 
rameters moves the gapless points in k-space, without 
removing them. 

In particular, when A„ = 0, the system is T-symmetric. 
The third condition in Eqs. (|B7j) then gives the wavenum- 
ber k to be one of the X r points; the other conditions 
lead us to the phase diagram in Ref. [l2| (also in Fig. [5Ja)) 
easily. 
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